Nonlinear resonance absorption in laser-cluster interaction 
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Rare gas or metal clusters are known to absorb laser energy very efficiently. Upon cluster ex- 
pansion the Mie plasma frequency may become equal to the laser frequency. This linear resonance 
has been well studied both experimentally and theoretically employing pump probe schemes. In 
this work we focus on the few-cycle regime or the early stage of the cluster dynamics where linear 
resonance is not met but nevertheless efficient absorption of laser energy persists. By retrieving 
time-dependent oscillator frequencies from particle-in-cell simulation results, we show that nonlin- 
ear resonance is the dominant mechanism behind outer ionization and energy absorption in near 
infrared laser-driven clusters. 

PACS numbers: 36.40.Gk, 52.25.Os, 52.50.Jm 



The construction of laser-based table-top sources of en- 
ergetic electrons, ions, and photons requires the efficient 
coupling of the incident (usually 800-1064 nm) laser light 
with matter. Clusters turned out to be very efficient ab- 
sorbers of laser light since they combine the advantages 
of gas and solid targets, namely transparency and high 
charge density. Almost 100% absorption of the laser light 
was observed in experiments with rare gas clusters [l|. 

The interaction scenario for laser wavelengths > 
800 nm, on which many researchers in the field agree 
upon, is as follows: after removal of the first electrons 
from their "parent" ions (inner ionization) and the clus- 
ter as a whole (outer ionization) the cluster charges up. 
The total electric field (i.e., laser plus space charge field) 
inside the cluster leads to inner ionization up to high 
charge states not possible with the laser field alone (ion- 
ization ignition 0,13 )• However, the restoring force of the 
ions counteracts outer ionization so that ionization igni- 
tion stops at some point. Moreover, the cluster expands 
due to Coulomb explosion and thermal pressure, thus 
lowering the electric field due to the ions. The latter de- 
termines the dominant eigenfrequency of the cluster, i.e., 
the Mie frequency couie(t) — k> p (t)/v3 = \/4np(t)/3 — 
y/Ni(t)Z(t)/R 3 {t) with p the charge density, N { the num- 
ber of ions of (mean) charge state Z, and R the cluster 
radius (atomic units are used unless noted otherwise). 
For laser wavelengths > 500 nm, soon after the removal 
of the first electrons WMie exceeds the laser frequency u>\. 
Hence linear resonance u>Mie(t) = w i occurs not before 
the cluster has sufficiently expanded (typically after a few 
hundred femtoseconds). At linear resonance the electric 
field inside the cluster is enhanced instead of shielded |4j 
so that even higher charge states can be created and even 
more energy can be absorbed from the laser. 

Absorption of laser energy is only possible through res- 
onances (linear or nonlinear) or non-adiabaticities (col- 
lisions). Commonly used phrases to explain absorption, 
like the tautological "laser dephasing heating" or the ill- 
defined "collisions with the cluster boundary," are correct 
but meaningless since dephasing is, according to Poynt- 



ing's theorem, a prerequisite for absorption while electron 
collisions with the "cluster boundary" occur all the time, 
independent of whether absorption is efficient or not. 

The importance of the linear resonance has been 
demonstrated both in pump pro be experiments and sim- 
ulations HH00liEalIJE3- The emission of 

third and higher harmonics of the incident laser light 
has been observed in computer simulations (l3|. Colli- 
sional ionization and absorption are of minor importance 
at wavelengths ~ 800 nm or greater |14l Il5l [l(| whereas 
it is the dominant absorption mechanism at short wave- 
lengths [lfll fl?l flij l not studied in the present Letter. 

One of the crucial points in the above mentioned sce- 
nario is the mechanism of outer ionization, which goes 
hand in hand with absorption |l9| since the laser en- 
ergy is transiently stored in the freed, energetic electrons. 
These electrons leave net positive charge behind, which 
finally Coulomb explodes. The latter converts electron 
energy into ion energy, which explains why experimental- 
ists typically measure MeV ions but only keV electrons 
[20L EH . In order to separate outer ionization and 
absorption due to linear resonance from other mecha- 
nisms we consider only the first few tens of femtoseconds 
of the laser-cluster interaction where ion motion is neg- 
ligible and linear resonance is therefore ruled out. Pre- 
vious work pointed out alr eady the possible importance 
of nonlinear resonance I23L |2 4l| or, equivalently, Landau 
damping in finite systems |25j. So far these mechanisms 
were clearly observed only in simplified model systems. 
The interpretation of results from molecular dynamics 
or particle-in-cell simulations, on the other hand, is of- 
ten hampered by the complex dynamics of the individual 
particles that makes the clear distinction of absorption 
mechanisms difficult. In this Letter we bridge this gap by 
analyzing our particle-in-cell results in terms of nonlinear 
oscillators. This enables us to prove that essentially all 
electrons contributing to outer ionization pass through 
the nonlinear resonance, so that the latter is unequivo- 
cally identified as the collisionless absorption mechanism 
if linear resonance is impossible. 
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FIG. 1: Laser absorption vs. driver amplitude in the 
rigid sphere-model @ for tt>Mie/wi = 2.7, -E(t) = 
— Ed sin 2 (r/2n) cos(r), n = 10. The nonlinear resonance is 
passed once the threshold driver amplitude ~ 3.0 is reached. 



In general, the eigenfrequency of a particle in a given 
potential depends on the excursion amplitude (or the en- 
ergy) of the particle, u> = w[f]. In a (laser-) driven system 
the excursion amplitude is time-dependent so that it may 
dynamically pass through the nonlinear resonance (NLR) 



w[r(t)] 



(1) 



Refs. [2J, |26j discuss at length the significant features of 
laser-driven systems undergoing NLR. Here we restrict 
ourselves to the prerequisites necessary to understand the 
analysis of our particlc-in-ccll results below. 

In the rigid sphere-model (RSM) (see, e.g., 0111 HI) 
of a cluster both electrons and ions are modelled by ho- 
mogeneously charged spheres which, in a linearly polar- 
ized laser field, oscillate along z around their common 
center of mass. In the case of equal charge density and 
radii R, the equation of motion for the center of the elec- 
tron sphere can be written in dimensionlcss entities as 



d r 



sgn(z) 



9r 2 
r 



r 
32 



Here, r = t — u>\t, and the upper line applies for 

< r < 2, the lower line (Coulomb force) for r > 2. As 
was observed in Ref. 0] the absorption of laser energy in 
the RSM is characterized by a threshold driver strength 
below which absorption is negligible (harmonic regime) 
and above which absorption is almost constant. Figure 
H shows this threshold behavior for WMic/^i = 2.7 and 
a n — 10-cycle sin 2 -pulse E(t) = — Eq sin 2 (r/2n) cos(r). 
Equation J2J can be formally rewritten as d 2 r/dr 2 + 



FIG. 2: (color online). Typical behavior of (o> e ff (t)/o>i) 
(drawn red) vs. time above the threshold driver strength in 



Fig. (actual value was Eq/Ru 



3.3). Excursion z/R 



(black) and energy of the electron sphere E to t/R 2 iof (green) 
are included in the plot. Outer ionization (i.e., E to t/R 2 u 2 > 
0) and occurrence of nonlinear resonance (u c ff (t)/lo\) 2 = 1 
(dashed-dotted line) always coincide (vertical line). 



(a; eff [f(r)]/a; 1 ) 2 sgn(z)r = E(t)/(R^) with 
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the instantaneous, scaled effective frequency lj c r(t)/uj\, 
which passes through unity at the NLR. Figure 
shows a typical example for the temporal behavior 
of (tJoff (t)/oj\) 2 above the threshold driver strength in 
Fig. Since (u> Mic /u)\) 2 = 7.29, (uj c g (t)/uj\) 2 starts at 
this value (dashed line in Fig. [21 and drops with increas- 
ing driver strength. It passes through unity at the time 
indicated by the vertical line, and it is exactly at that 
time where the electron sphere is set free, as it is clearly 
visible from the energy of the electron sphere, which 
passes through zero, and the excursion. We have checked 
numerically that NLR occurrs for all driver strengths 
above the threshold whereas the resonance is never met 
below threshold. In rare cases above the threshold, the 
electron sphere returns to the origin and recombines. In 
general, however, NLR results in a detached electron 
sphere at the end of the laser pulse. 

Let us now turn to the particle-in-cell (PIC) [2^ re- 
sults. We consider pre-ionized clusters of fixed radius 
R = 3.2 nm (e.g., Xe N , with Ni ~ 1600) but of vari- 
ous charge densities (i.e., different degree of inner ion- 
ization). The ratio of charge density to critical den- 
sity p/p c = 3 UJ Mic/ UJ \ var ies from 20 to 100. The clus- 
ters are exposed to 8-cycle sin 2 -pulses of wavelength 
A = 1056 nm. Since ion motion does not play an impor- 
tant role during the simulation time, the ions are fixed, 
which ensures a well defined, constant Mie frequency 
CiJMic- Electron-ion collisions are neglected in our PIC 
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FIG. 3: (color online). PIC results for a Xei6oo cluster (R = 
3.2 nm). Total absorbed energy per electron in units of U p vs. 
laser intensity for charge densities between 20 and 100 times 
the critical density, corresponding to average charge states 
between 1.6 and 8. The prediction of the RSM for p/p c = 40 
is included in the plot (dashed). 



treatment so that absorption of laser energy can only pro- 
ceed through other non-adiabaticities (like collisions with 
the "cluster boundary" or the cluster as a whole) or reso- 
nances. Figure 13 shows the absorbed energy per electron 
in units of the ponderomotive potential U p = Eq/Auj^, 
i.e., the time-averaged quiver energy of a free electron in 
the laser field. One sees that the absorbed energy per 
electron is always on the order of U p . However, the ab- 
sorbed energy is nonlinear in U p and displays a maximum 
before it drops because outer ionization saturates. The 
maxima are located close to the threshold intensities pre- 
dicted by the RSM (see RSM-result for p/p c = 40). With 
increasing charge density the maxima of the absorbed en- 
ergy (divided by U p ) move towards higher laser intensities 
while the absorbed energy per electron decreases. 

The motion of the PIC particles can be analyzed in 
the same way as it was done with the motion of the elec- 
tron sphere in the RSM above. A PIC particle has the 
same charge to mass ratio as a "real" electron, that is, 
e/m = —1 in atomic units. Each PIC particle moves 
under the influence of the external laser field and the 
space charge field E sc — — V ( f>(r,f) due to the potential 
$(r, t) that is created by all charges (mapped to a numer- 
ical grid). Hence the equation of motion of the ith PIC 
particle is + E sc (ri, t) = —E(t). The equation for the 
effective, time-dependent oscillator frequency analogous 
to then reads 



eff.iW — - 



[E(t) 



Esc(n,t) ■ Ti 



(4) 



Clearly, E sc (ri,t) depends on the position of all other 
particles ^ i as well. The PIC simulation starts with 
the neutral cluster configuration, i.e., the positive charge 



by the ionic background is neutralized by the negative 
charge of the PIC particles, so that E sc (ri,0) = 0. 
Hence, a PIC electron "sees" initially an effective fre- 
quency u> e e,i(0) = 0. The laser field disturbs the charge 
equilibrium and w 2 ffi (i) becomes different from zero. 
lo^ s . (t) may be even negative in regions of accumu- 
lated electron density (repulsive potential). As the clus- 
ter charges up, (w e ff/wi) 2 quickly increases beyond unity 
(where the RSM starts in the first place) . In the PIC sim- 
ulation each electron experiences its own u> e s,i, depending 
on the position. Hence, a sharp threshold in absorption 
as in Fig.^does not exist. Besides this, the starting from 
w e ff,i(0) = 0, the possibility of negative w 2 ff and the 
three-dimensionality are the main differences to the RSM 
analysis above. 

By following the dynamics of the electrons in the ef- 
fective frequency vs. energy-plane one can identify the 
main pathway to outer ionization and efficient absorp- 
tion. In Fig. 0Ji-d the scaled effective frequencies squared 
(w e ff/uJi) 2 of the individual PIC electrons are plotted 
vs. their energies E totyi (t) — r?(i)/2 — $(rj,i) they 
would have if the driver is switched off instantaneously 
at t — 2.5, 3, 3.5, and 4 laser cycles, respectively. We 
define the time when, for a particular electron, E to t be- 
comes > as the ionization time of that electron. The 
laser intensity is 2.5 x 10 16 Wcm~ 2 , and the preionized 
cluster is 40 times overcritical so that (wMio/^i) 2 =40/3. 
As is clearly visible in Fig. 0] each electron reaches posi- 
tive energy close to the point (w 2 ff /w 2 , E tot /U p ) = (1, 0). 
The radial position of each electron is color coded, indi- 
cating that outer ionization occurs at radii around 2R. 
Data points at radii > 3i? (orange and yellow colors in 
Fig. Et,b) with positive but very small E tot and w 2 ff ~ 
represent low energetic electrons removed earlier during 
the pulse. Electrons with positive energy but small radii 
[visible in (a) and (b)] are those driven back to the cluster 
by the laser field. For the electrons inside the cluster po- 
tential (negative energies and radii < R, color coded blue 
and black) (uj c g/uji) 2 spreads over a wide range, starting 
from the maximum value (wMie/^i) 2 down to negative 
values due to the repulsive force exerted by the com- 
pressed electron cloud. Note that negative values occur 
mainly at early times where most of the electrons are still 
inside the cluster. The occurrence of NLR is less clear 
for the few electrons leaving the cluster even earlier than 
t = 2.5 laser cycles. As mentioned above, these electrons 
move in a shallow effective potential with (w c ff/wi) 2 < 1 
when they leave the cluster with ease and with rather 
low kinetic energy because the laser intensity is still low 
at the time of their emission. With increasing degree of 
outer ionization, the resonance point (1,0) is met more 
and more precisely [i.e., less vertical spread around (1,0)] 
since the effective frequency is more and more dominated 
by the uniform, static ion background. The electrons 
pass very quickly through the resonance (typically less 
than a quarter of a laser cycle) because the energy gain 
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FIG. 4: (color online). Snapshot of PIC 
quency vs. energy-plane at times (a) t 
(c) £ = 3.5, and (d) t = 4.0 laser cycles, 
was 2.5 x 10 16 Wcm _2 and (aj M ic/^i) 2 = 
eters as in Fig. |3] The radial positions 
color-coded. Electrons become free upon 
ear resonance, i.e., (u> 2 ff /w 2 , E to t/U p ) = 
the radial position is r ~ 2. 



electrons in the fre- 
= 2.5, (b) t = 3.0, 

The laser intensity 
40/3. Other param- 

(in units of 7?) are 
crossing the nonlin- 
(1,0). At that time 



necessarily drives the electrons out of resonance again. 

One may object that, since the denominator in Q 
necessarily increases while the numerator decreases for 
an electron on its way out of the cluster potential, the 
passage through a point (iv^ s /ujf, E tot /U p ) = (x,0) with 
x some value < (wMie/wi) 2 is rather the consequence of 
outer ionization than the mechanism behind it. How- 
ever, NLR only occurs at x = 1, and the results in Fig.0] 



show only little spreading along (w e ff/wi) at E tl 



0. 



Moreover, the fact that both the single electron ener- 
gies become positive and the radii exceed ~ 2R when 
(w e ff/u;i) 2 = 1 indicates that NLR is indeed the responsi- 
ble mechanism behind outer ionization accompanied by 
efficient absorption of laser energy. One may further ob- 
ject that such a short resonance is more like a collision 
(e.g., with the cluster boundary) than a resonance. How- 
ever, note that collisions with the cluster boundary occur 
many times but efficient absorption and outer ionization 
only takes place when the resonance condition is fulfilled. 

In summary, we have shown that cluster electrons con- 
tributing to efficient absorption and outer ionization in 
near infrared laser fields undergo nonlinear resonance, 
meaning that the instantaneous frequency of their mo- 
tion in a time-dependent, anharmonic, effective poten- 
tial meets the laser frequency. Nonlinear resonance is 



the only possible absorption mechanism if the laser pulse 
is too short for the linear resonance to occur (or dur- 
ing the early cluster dynamics in longer pulses) and if 
electron-ion collisions (inverse bremsstrahlung) are neg- 
ligible, as it is the case at near infrared or greater wave- 
lengths. In order to prove the occurrence of nonlinear 
resonance we introduced a method to analyze the re- 
sults obtained from particle-in-cell simulations, namely 
the mapping of the system of electrons and ions that 
interact through their mean field onto a system of non- 
linear oscillators whose time-dependent frequencies un- 
equivocally revealed the coincidence of electron removal 
and nonlinear resonance. 

Fruitful discussions with P. Mulser are gratefully ac- 
knowledged. This work was supported by the Deutsche 
Forschungsgemeinschaft . 
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